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ABSTRACT 



Aims. We use the super Earth 55 Cnc e as a case study to address an observable effect of tidal heating. This transiting short-period 
planet belongs to a compact multiple system with massive planets. We investigate whether planet-planet interactions can force the 
eccentricity of this planet to a level affecting the eclipse depth observed with Spitzer. 

Methods. Using the constant time lag tidal model, we first calculate the observed planet flux as a function of albedo and eccentricity, 
for different tidal dissipation constants and for two extreme cases: a planet with no heat redistribution and a planet with full heat redis- 
tribution. We derive the values of albedo and eccentricity that match the observed transit depth. We then perform N-body simulations 
of the planetary system including tides and General Relativity to follow the evolution of the eccentricity of planet e. We compare the 
range of eccentricities given by the simulations with the eccentricities required to alter the eclipse depth. 

Results. Using our nominal value for the dissipation constant and the most recent estimates of the orbital elements and masses of the 
55 Cnc planets, we find that the eccentricity of planet e can be large enough to contribute at a measurable level to the thermal emission 
measured with Spitzer. This affects the constraints on the albedo of the planet, which can be as high as 0.9 (instead of 0.55 when 
ignoring tidal heating). We also derive a maximum value for the eccentricity of planet e directly from the eclipse depth: e < 0.015 
assuming Earth's dissipation constant. 

Conclusions. Transiting exoplanets in multiple planet systems - like 55 Cancri - are exceptional targets for testing tidal models be- 
cause their tidal luminosity may be observable. Future multi- wavelengths observations of eclipse depth and phase curves (for instance 
with EChO and JWST) should allow us to better resolve the temperature map of these planets and break the degeneracy between 
albedo and tidal heating that remains for single band observations. In addition, an accurate determination of the eccentricity will make 
it possible to constrain the dissipation rate of the planet and to probe its internal structure. 



Key words. Planets and satellites: fundamental parameters - 
interactions - 



1. Introduction 



The bright star 55 Cancri hosts at least five planets (Fis cher etal. 
2008 ) including t he close-in 8M @ , 2. 177? @ planet e that is on a 



transi ting orbit (D emory et al) |2011[ |Winn et al.| |2011| |Gillon 
|et al.||2012| ). Table [T] shows the parameters of the star and planet 
e. Planets 55 Cnc e and Kepler 10 b are the only low-mass plan- 
ets (M < 10M e ) whose seco ndary eclipse has been observed, 
respectively by Demory et al. ( |2012| hereafter D12) and Batalha 
et al. ( |2011| ). In the case of Kepler 10 b, the planetary signal re- 
sponsible for the eclipse depth is dominated by reflected light 
due to the spectral window of Kepler (400 - 900 nm). For this 
planet, a high albedo (A > 0.5) is required to fit the observation 
( [Batalha et al.||201l||Rouan et al.||2011| ). The secondary eclipse 



Planets and satellites: dynamical evolution and stability - Planet- star 



put tight constraints on the radius of 55 Cnc e by combining 
primary transit data from both MOST and Spitzer. The parame- 
ters of the bright host star 55 Cnc have been accurately charac- 
terized using various techniques (including interferometric mea- 
surements) by |von Braun et al.| ( |201 1\ . 



of 55 Cnc e has been measured by D12 using the 4-5 jum IRAC2 
channel of Spitzer. At these wavelengths, the thermal emission 
dominates the planetary flux. The observation implies an upper 
limit on the albedo (A < 0.55) as well as an inefficient heat trans- 
port towards the night side of the planet. Even with a null albedo, 
a planet with a uniform temperature would have an insufficient 
dayside emission to produce the observed drop of flux during the 
secondary eclipse (D12). The required albedo depends of course 
strongly on the radius of the planet but Gillon et al. (2012 ) have 



Planets Kepler 10b and 55 Cnc e share a lot of similarities 
but their surface albedo derived from secondary eclipse seem 
to differ. The composition of the two planets can of course be 
different. In particular, 55 Cnc e has a surprisingly low den- 
sity that implies a significant volatile content (and therefore a 
gaseous envelope), as suggested by D 12, or a carbon-rich com- 
position (ceramics being less dense than silicates) as suggested 
by Madhusudh an et aL] ( |2012| ). In the present study, we inves- 
tigate a scenario that could change the albedo determination, in 
which the thermal emission of 55 Cnc e is enhanced by the tidal 
dissipation resulting from a low eccentricity maintained by the 
other pla nets of the system. With a method similar to the one 
Barnes et al. (2010) used to infer the eccentricity of Corot-7 b 
and its tidal heating, we used N-body simulations that include 
tidal interactions and General Relativity to study the evolution 
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8.26 M e 


2.11 Re 


0.01560 









Table 1. Stellar properties and minimum mass, radius, and semi-major axis of planet 55 Cnc e. The stellar mass, radius, effective 
temperature, gravity and metallicity come from |von Braun et aL](|201 1|) and its rotation period comes from |Fischer et al.|(|2008|). Th e 
values of M sin / and a come from table 10 of Dawson & Fab ry cky (2010). The radius of 55 Cnc e comes from Gillon et al. 12012). 



of the system and of 55 Cnc e in particular. We find that the orbit 
of 55 Cnc e cannot be fully circularized and keeps a low ec- 
centricity subjected to important variations. Using the planets' 
orbital elements derived from radial velocity data, we find that 
the eccentricity of planet e can reach values close to 0.01, which 
is below the upper limit from observations (e=0.06, D12) but 
enough for dissipation to produce an internal heat flux compara- 
ble to the stellar heating. In this case, the tidal dissipation affects 
the observed secondary eclipse depth, which can be compatible 
with a high albedo. 

In a first section we study the influence of tidal dissipation on the 
thermal emission of the planet and derive the minimum eccen- 
tricities required to affect the thermal phase curve and secondary 
eclipse of 55 Cnc e, for two extreme cases: a planet with no heat 
redistribution - no atmosphere - and a planet with a full redistri- 
bution of the stellar heating. In a second section, we simulate the 
dynamical evolution of the 55 Cnc system in order to constraint 
the range of possible eccentricities for planet e. 

2. Effect of eccentricity on the secondary transit 
depth. 

2.1. Modeling the planet flux 

In order to bracket the range of possible situations we consider 
two different cases: a synchronously-rotating planet with no at- 
mosphere and a planet with a uniform temperature due to a fully 
efficient redistribution of the stellar heating by an atmosphere. 
The bolometric stellar flux 0* and the spectral flux density 0* >/t 
received at a distance d are given by given by 



b*(d) = o-T 



if 



m 



where cr is the Stefan-Boltzmann constant, T e ff and R+ are the 
effective temperature and radius of the star, and 



b* A {d) = F+ A 



(#• 



where F*^ is the spectral flux density at the surface of the 
star and is computed with the stellar atmosphere code Phoenix 
( |Allard et all |2012) for the T eff , gravity (log(g)) and metallic- 
ity of 55 Cnc (see Table [T]). This synthetic spectrum agrees well 
with available observations (I. Crossfield, private communica- 
tion), which are numerous for such a bright star (e.g. IRAS, 
IRTF, Akari, Spitzer, Wise). 

The planet spectrum is a combination of its thermal emis- 
sion and the reflected light. To derive the thermal component we 
first need to calculate the surface temperature map, which de- 
pends on both the absorbed stellar flux and the internal heat flux. 
Note that we consider an eccentricity that remains small enough 
so that its effect on the irradiation remains negligible. To calcu- 
late the stellar heating, we can thus assume a circular orbit and a 



synchronous rotation. In those circumstances, the irradiation at 
a given point of the planet does not vary with time. We checked 
with a more detailed thermophysical model that accounting for 
the small librations of the planet and the thermal inertia of the 
surface does not change our results. 

If we assume no redistribution of the incident flux (no atmo- 
sphere), then the stellar flux absorbed by the planetary surface is 
a function of the zenith angle 6 and the Bond albedo A B and is 
given by 

0abs = 0*(a)cos0(l -A B ), 

a being the semi-major axis of the orbit. As we set the ori- 
gin of the geographic coordinates at the substellar point, the 
zenith angle is given by cos = cos(latitude) cos(longitude). We 
also include an internal heat flow due to tidal dissipation. Other 
sources of internal heat like radioactivity or remnant of accre- 
tion/differentiation can dominate when dissipation is negligible 
but in this case the resulting surface warming would be too low 
to affect the observed emission. For instance, the mean heat flow 
of the Earth corresponds to a surface temperature of only 30 K. 
Significantly higher values would imply a very young planet, a 
case we do not consider here given von Braun et al. ( 201 1 )'s esti- 
mation of the age of 55 Cnc: 10.2 Gyr. Locally, the flux emitted 
by the surface is equal to the sum of the absorbed stellar flux 
and the internal heat flow t id es resulting from tidal dissipation, 
which we assume uniform over the whole planetary surface. On 
the dayside, the surface temperature is thus given by 



?abs 



(fl) + , 



'tides 



while on the night side we have 

^night = (0tides/0") 5 - 

In the case of a full redistribution of the absorbed stellar en- 
ergy the temperature is uniform over the whole planet and sim- 
ply given by 

where T Qq is the equilibrium temperature of the planet: 
(</>*(a)(l-A B )\* 



^eq — 



Act 



The surface of the planet is divided into a longitude-latitude grid. 
Each cell j of the grid has a surface temperature I) and an area 
Sj. The flux spectral density received by a distant observer at a 
distance d is given by 



(f>P,A(d) = Y J I J — 



S j cos a j 



d 2 
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where aj is the angle between the normal to the cell and the 
direction toward the observer and Ij is the specific intensity of 
the cell given by 

_ e A B A (Tj) + A x <p* A (a) cos 0,- 



where B A is the Planck function, e A is the surface emissivity, and 
A A is the surface spectral albedo. In this study, we use e A = 1 and 
A A = A B = A. In practice, cos aj is calculated as the dot product 
of the directions vectors attached to the center of the planet and 
pointing towards the observer and the center of the j cell. Only 
locations visible to the observer (cos aj > 0) contribute to 0/^. 

In the previous formula, the star is considered as a point 
source but this approximation can become inadequate in the case 
of extremely close planets. Seen from 55 Cnc e, the angular ra- 
dius of the stellar disk is 15°. This has two main consequences. 
First, it affects the distribution of the stellar flux over the plan- 
etary surface and thus the temperature map. The terminator that 
delimits the dayside and the nightsi de is re placed by a penumbra 
region, as described by Leger et al. ( [201 1] ). This effect, however, 
has a negligible impact on the observed flux even in the case of 
55 Cnc e because the affected regions of the planet emit/reflect 
an insignificant fraction of the total spatially-unresolved plane- 
tary flux. We simulated numerically an extended star by spread- 
ing the stellar luminosity over a number of fainter punctual 
sources distributed over the stellar disk. We found that, for phase 
angles smaller than 120°, the planetary flux calculated with a 
point source star and with an extended star differ by less than 1 %. 
There is a second, more significant, effect: phase angles smaller 
than 15° cannot be observed as they are eclipsed by the star. The 
out-of- transit photometric reference used by D12 actually cov- 
ers a broad range of phase angles: 68 - 15° before the eclipse 
and 15-43° after, while the eclipse itself lasts about one tenth 
of the orbit. With our model, we can calculate the planetary flux 
from an airless planet for any phase angle (i.e. sub-observer lon- 
gitude). Therefore, we can either produce the light curve for the 
whole range of observed phase angles (as in Fig. [T] to [2]) or for 
a phase angle of 35°, which is the averaged value of the out-of- 
transit observations by D12 (as in Fig[3]to|4]). 

To compare the computed planetary flux with the measured 
transit depth of 131 ±28 ppm derived by D12, we calculate the 
planet to star contrast ratio R as follows: 



R = 



j AA (Pp,a(D)w a dA 
f M (f>*, A (D)w A dA 



where D = 12.34 pc is the distance of the 55 Cnc system, w A is 
the spectral response of the IRAC2 detector (as given in Fig. 2.4 
of the IRAC Instrument Handbook), and AA is the IRAC2 band 
(3.9- 5.1 //m). 

2.1.1. Tidal heating 

The average surface heat flux t i des is given by E ti( i e sl47iR 2 p where 
Etides i s me rate of tidal dissipation. We use here the constant 
time lag model to compute E tides- For a planet of mass M p in 
pseudo-synchronous rotation orbiting a star of mass M* at a 
semi-major axis a and with a null obliquity (formula 13 from 
|Leconte et al.||201Q| ), the rate of tidal dissipation is given by: 



Etide ~ 2 
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Fig. 1. The effect of albedo and eccentricity on the light curve. 
This graph shows the planet/star contrast ratio as a function of 
the orbital phase for different values of (e, A). All these models 
give the same flux (131 ppm) for a phase angle of 35°, which is 
the average value for the out-of-eclipse Spitzer observations by 
D12. Calculations are done with the nomina l value cr p and th e 
orbital parameters and planetary mass from E ndl et al.| ( |2Q12| ). 
The short-dashed and full lines are for a uniform temperature 
(perfect redistribution of heat) while there is no redistribution in 
the three other cases. 
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where the dissipation timescale T p is defined as 



1 



M n 



1 



9M*(M p + M*)fli°o> 



(2) 



and depends on the mass of the planet M p , its dissipation factor 
(T p (as defined by Hansen (2010)) and of the mass of the star M*. 
Nal(e), Na2(e) and Ll(e) are eccentricity-dependent factors: 



Nal(e) 

Na2(e) 

0{e) 



1 + 31/2e 2 + 255/8e 4 + 185/16e 6 + 85/64<g 8 

(1 - e 2 ) 15 ' 1 
l + 15/2e 2 + 45/8<? 4 + 5/16e 6 
(1 - e 2 ) 6 ' 

1 + 3e 2 + 3/8e 4 
(1 - e 2 ) 9 ' 2 ' 
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Fig. 3. Contraints on albedo and eccentricity from the secondary eclipse depth, in the case of an airless planets. The light green 
shaded areas are in agreement with the 131 ±28 ppm depth measured with IRAC2/Spitzer by D12. The right hand side plot is for 
an Earth-like dissipation constant cr p ; the left and middle plots are for 10" 3 x cr p and 10" 2 x cr p , respectively. The dotted (resp. 
dashed) lines indicate where reflection (resp. tidal dissipation) accounts for 10 and 50% of the planetary flux (reflection + thermal 
re-emission + tides = 100%). The long dashed pink vertical line represents the observational constraint on the eccentricity given by 
D12. Horiz ontal grey bars below the graphs indicate the range of eccentricities found with the dy na mical simulations described in 

1( [2QT2l > and |Dawson & Fabrycky| ([2010), 



3.2 The upper and lower bars correspond to initial orbital elements from Endl et al 



le minimum when larger than 5 X 10 



section 

respectively. The maximum and mean values of the eccentricity are indicated, as well as th 

The blue areas are in agreement with the transit depth, the condition that tides accounts for more than 10% of the planetary flux and 

compatible with an eccentricity within the range found with the dynamical simulations. For 10~ 3 xcr p and 10~ 2 X cr p the agreement 

f [pT..,. _ o, n,u ,,., ,™„y, ...u:,.. t. , I. , .u ,, -^ + >... r, „ ,, „n, 



is met for simulations of Dawson & Fabry cky (2010) while for 1 x cr p the agreement is met for Endl et al. (2012 ). 



The difference in the rate of tidal dissipation assuming syn- 
chronization or pseudo- synchronization for planet 55 Cnc e is 
3.5% for e = 0.06 (the upper limit of eccentricity from D12) and 
only 0.1% for e = 0.01. 

The characteristics of planet e can be seen in Table [T] We 
used for this study the lat est est imation for the radius of 55 C nc e 
from|Gillon et al | { |2012) . Both pawson & FabryckylpOTO] ) and 
Endl et al. (2012) give comparable masses for 55 Cnc e. 

As in Bolmont e t al.| ( |2011| |2012| ) in order to be consistent 
with the Ea rth dissipation i nferred from D245 astrometric data 
(Neron de Surgy & Laskar 1997) ), our nominal value of cr p is 
1.877 x 10" M g _i cm~ 2 s _1 . This value will be referred thereafter 
as 55 Cnc e nominal dissipation factor. 

The uncertainty on the dissipation factor is, however, consid- 
erable. The composition of the planet and its internal structure 
are unknown while they affect, and are affected by, the dissipa- 
tion. For instance, if we assume that a melted planet dissipates 
less than a solid one, an initially solid planet can start to melt 
due to tidal heating, which in turn lowers its dissipation and the 
internal heating. For a given eccentricity, this feedback may re- 
sult in a steady state, in which the internal structure and mineral 
phases are consistent with the tidal heating. A self-consistent 
treatment would thus include a feedback between the dissipation 
rate and the dissipation factor, but that is beyond the scope of this 
stu dy. (Such modeling has be en done by |Henning et al.| ( |2009| ) 
and Behou nkova et al.| ( |2011| ).) Nonetheless, there is no reason 
that 55 Cnc e's dissipation factor should be the same as Earth's. 
To take into account this uncertainty we performed simulations 
over 7 decades of dissipation factors from 10" 5 to 10 2 x cr p . 



Figure [5] shows the contours of the average surface heat flux 
0tides> as a function of eccentricity and dissipation factor. To af- 
fect the observation, t id es must be extremely high and compara- 
ble to the insolation flux. Such values exceed by orders of mag- 
nitude the tidal heat flux measured at the surface of Io, which 
is the largest in the Solar System. The highest values on Fig. [5] 
certainly imply a melted surface and interior, as well as an ex- 
tremely efficient vertical heat transfer, which may not be con- 
sistent with the assumed strong dissipation, as discussed above. 
Note, however, that a melted surface does not contradict the air- 
less assumption as melted refractory materials have very low va- 
por pressures: Lege r et aL] ( |2011| ) calculated that the equilibrium 
pressure above a substellar magma ocean on Corot 7 b would 
not exceed a few Pa, with a negligible effect on the opacity and 
horizontal heat transport. 

2.2. Model vs observations 

Besides parameters that are constrained by observations other 
than secondary eclipse (M+,Mp,Rp,d), t ides depends on 
the eccentricity and on the internal dissipation constant cr p , 
while the albedo A determines the ratio between reflected and 
absorbed stellar fluxes. At low eccentricities yielding negligible 
dissipation, the planetary flux depends only on the albedo (and 
on the heat redistribution assumption). For higher eccentricities, 
the resulting dissipation can increase the temperature and the 
flux depends on two additional parameters: e and cr p . 
Figures [3] and [4] show contours of the eclipse depth in the 
IRAC2/Spitzer band as a function of eccentricity and albedo, 
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Fig. 4. Same as in Fig.[3Jbut for an isothermal planetary surface. 
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Fig. 5. Map representing contours of the average surface heat 
flux t ides> as a function of eccentricity and planetary dissipation 
factor. From left to right the full lines correspond to a flux of 
10" 6 , 10" 4 , 10" 2 1, 10 2 , 10 4 , 10 6 , 10 8 , 10 10 W.m" 2 . The long 
dashed lines correspond to a temperature T = (fades /v)* of 
10, 10 2 , 10 3 a nd 10 4 K. The obse rved value of t id es on Io, 
2.4-4.8 W/m 2 ( Spencer et~aTl|2000| ), is ind icated in black for its 
estimated eccentricity value of 4.3 x 10" 3 ( Murray & Dermott 
[T999] ). 



and for three values of the dissipation constant: 10 3 x cr D 



10" 2 x cr p and 1 x a_p. Figure p\ is for a planet with no at- 
mosphere while Fig. [4] is for an isothermal photosphere (full 
redistribution of the stellar heating). The shaded area gives 
the A — e domain compatible with the observation by D12 
(131 ± 28 ppm). The eclipse depth is calculated for a phase 
angle of 35° corresponding to the average phase angle during 
the out-of-eclipse observations by D12. Isothermal cases cannot 
match the observation even for A = (as noted by D12) because 



the dayside would be too cold and would produce an eclipse 
that is too shallow unless there is an internal source of heat. 
This isothermal case may not, however, be realistic for such 
a close-in planet, even in the presence of a thick atmosphere. 
Indeed, the radiative cooling timescale, which varies as T~ 3 , 
is likely to be shorter than the dynamical one due to the high 
equilibrium temperature of this planet (~ 1600 K, for an albedo 
of 0.5). Even a moderate day-night temperature contrasts would 
make the dayside emission significantly stronger than the 
nightside one. Note also that a phase modulation exists even in 
the isothermal case due to the reflected light (unless A = 0). 

The cases with no heat redistribution match the observation 
with either an albedo lower than 0.55 and no eccentricity (in 
agreement with D12) or with a higher albedo (up to A = 1) 
and an eccentric orbit. When increasing (resp. decreasing) the 
dissipation constant, lower (resp. higher) values of the eccen- 
tricity are required to modify the eclipse depth by a same fac- 
tor. In Fig. [3] and [4] the blue shaded area corresponds to cases 
for which the following three conditions are met: the resulting 
transit depth is within the observational range, the tides account 
for more than 10% of the planetary flux and the eccentricity is 
within the range of eccentricities found by the dynamical simu- 
lations of the system. For dissipations of 10" 3 xcr p and 10" 2 xcr p , 
the range of eccentricity comes from a simulation of the system 
of pawson & Fabry cky (2010) and for a dis sipation of 1 x cr p , 
it comes from a simulation of the system of Endl et al. (2012). 
For example, assuming that tides represent more than 10% of 
the planetary flux, the observations would be in agreement with 
a planet of albedo between 0.25 and 0.75 assuming a dissipation 
of 10" 3 x <j p . As we do not know the dissipation of the planet 
nor its precise eccentricity the observations would not allow to 
determine the albedo of the planet. However, D12 gave an obser- 
vational constraint on the eccentricity of 55 Cnc e, it should not 
be higher than 0.06. It allows to rule out some of the solutions 
found. For a dissipation of 10 ~ 3 x <t p , there would no longer be 
agreement with Dawson & Fabry cky (2010) because the eccen- 
tricity is too high. For a dissipation of 10~ z x cr p , there would be 



an agreement with pawson & Fabry cky (2010) for a planet of 
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albedo between 0.35 and 0.7. Fo r a dissipation of 1 x cr p 



would still be an agreement with Endl et al. ( 20 1 2| ) for a planet 
of albedo lower than 0.95. 

We define here the critical eccentricity as the eccentricity 
above which tides affect the transit depth by more than 10% 
and for which the transit depth is of 103 ppm. This critical 
eccentricity depends on the dissipation and is given in Table [5] 
for the range of dissipation factors considered here. In Figures 
[3] and [4] the critical eccentricity corresponds to the eccentricity 
of the intersection of the dashed line corresponding to "tides 
= 10%" and the upper line of the shaded area. The maximum 
eclipse depth that we find without dissipation is 134 ppm 
with no atmosphere and 98 ppm in the isothermal case (to be 
compared with the 131 ± 28 ppm found by D12). 

Whatever the redistribution efficiency, the albedo, and for a 
dissipation constant of 1 x cr p , eccentricities higher than 0.015 
produce too large a planet flux and can be ruled out. However, 
D12 constrained the eccentricity of 55 Cnc e to be less than 
0.06, which allows us to rule out some of the solutions we found. 

Fig. [3] and [4] also show the range of eccentricities found with 
N-body simulations that are described in the next section. These 
N-body simulations use two different sets of orbital elements 
and masses derived from radial velocit ies, one from |Dawson| 
|& Fabrycky] pOTO] ) and the other from |Endl et al.| (201 2) . In 
particular, the value for the orbital distance of planet e slightly 
differs between these two sets (see Tabl^3j. The effect on the 
transit depth calculation is, however, very small and would 
hardly be seen on these plots. Therefore, Fig. [3] and [4] only 
present the results obtained with a = 0.01544 AU, the value 
from |Endl etaTl { |20l2) . 

As shown in Fig.[T] orbital photometry (i.e. phase curve mea- 
surement) could break the degeneracy between the different sets 
of albedo and eccentricity. Cases with strong tidal heating and 
high albedo exhibit much flatter phase curves than cases with 
negligible tidal heating and low albedo. Distinguishing between 
these cases at 4.5 /mi would require a photometric precision of 
the order of 10 ppm, which is not achieved yet as seen in Fig. [2] 
The space telescope EChO ( [Tinetti et al.| |2012| ) should reach 
the required precision and could also observe at wavelengths as 
long as 11 yum and possibly 15 yum. At these wavelengths the 
planet/star contrast ratio would increase significantly (200 ppm 
at 1 1 /mi). Also by measuring the eclipse depth and/or the phase 
curve at different wavelengths simultaneously, EChO could get a 
much better constraint on the temperature map of the planet and 
thus on the albedo vs eccentricity. 

With no other constraints on the eccentricity, our results 
show that all values of the albedos between and 1 are possi- 
ble for 55 Cnc e if its dissipation factor is of 10" 2 x cr p . The goal 
of the next section is to estimate the maximum eccentricity of 
planet e using N-body simulations of the 55 Cnc system. 



3. The eccentricity of 55 Cnc e: damping vs forcing 

In a system with one planet, the eccentricity of such a close-in 
planet would be damped to zero very quickly, however 55 Cnc e 
is part of a compact five-planet system. Planet-planet interac- 
tions can maintain an eccentricity, which would be an equilib- 
rium between the tidal damping and the excitation by the other 
planets. 



there 3.1. Eccentricity evolution assuming a single planet system 



The model used for the study of the evolution of one planet 
around it s host star i s a a re-derivation of the eq uilibrium tide 
model of Hut (1981 ) as in Eggle ton et al.| ( [T998] ). We consider 
both the tide raised by the star on the planet and by the planet on 
the st ar. We apply the c onstant time lag model ( [Leconte et al. 



|2010| ) and use Hansen's ( |2010| ) estimation of the internal dissi- 
pation constant of giant exoplanets and their host stars for planet 
'b' and 55 Cnc. For consistency, we use for planet e the sa me dis - 
sipation constant as in the tidal heating model (see section 2.1.1). 
The secular evolution of the eccentricity of the planet is 
given by|Hansen|(|2010|): 



ld£?i ldei 

e dt \plan e dt \star 



2T X 



r 11 Hp 

[Nel(e)- — —Ne2(e) 



18 n 



(3) 



27V 



[Nem-^Ne2(e)\ 



where Q p is the planet rotation frequency, and n is the mean 
orbital angular frequency. The stellar parameters are obtained 
by switching the p and * indices. Nel(e) and Ne2(e) are given 
by: 



Ne\(e) = 
Ne2(e) = 



1 + 15/4e 2 + 15/8e 4 + 5/64e 6 

(1 - e 2 )^ 2 
l + 3/2e 2 + l/8<* 4 
(1 - e 2 ) 5 ' 



The integration of the tidal equations was performed using 
a fourth order Runge-Kutta integrator with an adaptive timestep 
routine (Pr ess etal.|| 1992] ). The precision of the calculations was 
chosen such that the final semi-major axis of each integrated sys- 
tem was robust to numerical error at a level of at most one part 
in 10 3 . 

The characteristics of the star and the inner planet of the 
55 Cnc system are listed in Table [T] 

We use for the s tar 55 Cnc the value of stellar dissipation 
derived by 



Hansen 



(2010): a 



4.992 x KT* g 



-1 -2 - 

cm s 



We compute the tidal evolution of 55 Cnc e with the dissipa- 
tion factor given in section 2.1.1 and we also compute sepa- 
rately the tidal evolution of 55 Cnc b. This planet is thought 
to be a gas giant, so Hansen's value for gas giants is used: 
cr b = 2.006 x 10" 60 g^cnrV 1 . The radius of 55 Cnc b is cal- 
culated assuming a Jupiter density. 

Because the planetary spin synchronization timescale is 
shor t compar ed to the other time scales considered here (Leconte 
|et al.| |2010[ jHeller et al!) |2011| ), each planet rotation period is 
fixed to the pseudo- synchronization value at every calculation 
timestep. Note that for the low eccentricities found for planet 
e, the difference in the dissipation rates calculated assuming 
pseudo-synchronization and sychronization is negligible. If the 
obliquity of planet e is determined by planetary tides, it reaches 
zero on very short timescales. The obliquity of planet 'b' should 
not evolve tidally on less than Gyr timescales. For simplicity, we 
assume both obliquities to be zero. 

If the obliquity of 55 Cnc e were non-zero it would affect our 
results by increasing the planet's dissipation rate and tidal heat 
flux. Neglecting the obliquity of 55 Cnc e therefore gives us a 
lower limit on the possible tidal flux. We calculated the tidal flux 
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for two different obliquities of planet e. Fig. [6] shows the ratio of 
the planetary tidal- to total flux in the Spitzer band for different 
obliquities, for dissipation factors and albedos. As expected, the 
lower the dissipation factor, the higher the eccentricity needed to 
reach the ratio of 10% of the planet tidal flux over its total flux. 
For a dissipation of 10" 3 x cr p , the eccentricity needed to reach 
10% is of about 0.1 for a planet with an albedo of and this 
eccentricity is of 0.04 if the albedo is of 1. In the IRAC2 Spitzer 
band, the thermal emission dominates over the reflected light so 
increasing the albedo has an effect of decreasing the non tidal 
part of the flux emitted by the planet. Thus, a smaller eccentricity 
is needed for the tidal flux to reach 10% of the planet total flux. 

The right panel of Fig. [6] shows a case for which planet e's 
obliquity of 2 degrees may have been excited by forcing from 
the other planets. Had this planet been alone, tides would have 
reduced this obliquity in very short timescales. The major differ- 
ence with the left panel is that for small eccentricities, the ratio 
^p,tides/^p,tot is non zero. The obliquity contributes to the tidal 
dissipation in a strong way. If the planetary dissipation is more 
than 10" 2 x cr p , the ratio is always bigger than 10% whatever the 
value of the planet's albedo. For eccentricities bigger than 0.1, 
the tidal dissipation due to the eccentricity becomes dominant 
again. 

The study of the effect of obliquities of the inner planets on 
the N-body evolution of the system will be conducted in the fu- 
ture. 



obi = deg 



For 55 Cnc e, r Pi 



4600 yrs, and the eccentricity damp- 



1.4 x 10 10 




Fig. 6. Ratio of 55 Cnc e tidal flux over its total flux as a func- 
tion of eccentricity for two different obliquities (0° and 2°). The 
purple curve corresponds to a dissipation factor of 10" 3 x cr p , the 
blue one to 10" 2 x cr p , the green one to 10" 1 x cr p and the orange 
one to 10 1 x cr p . The full line corresponds to a Bond albedo of 
the planet of 0, the small dashed line to an albedo of 0.5, the 
dashed dotted line to an albedo of 0.8 and the long dashed one 
to an albedo of 1 . The black horizontal dashed line represents a 
ratio of 10%. 



ecc \plan 

ing timescale due to the stellar tide is r ecc \ ~ 1.4 x 1U W yrs. 
The evolution of the eccentricity of 55 Cnc e will therefore be 
dominated by the planetary tide. Let us define the eccentricity 
damping timescale as r ecc = r ecc \ . In Fig.p7J the evolutionary 

tracks of 55 Cnc e are represented in purple aashed-dotted lines 
for a nominal dissipation factor cr p . If this planet was alone in the 
system, the eccentricity would be damped in about 10 x r ecc \ 

For 55 Cnc b, r ecc \ , and r ecc \ , are longer than the age 
of the system. The eccentricity of 55 Cnc b will remain at its 
present values on more than Gyrs timescales. In Fig. [7] the evolu- 
tionary tracks of 55 Cnc b are represented in blue dashed-dotted 
lines. If this planet was alone in the system, its semi-major axis 
and eccentricity would remain the same for more than one bil- 
lion years. 

However, to accurately describe the evolution of 55 Cnc e, 
it is crucial to take into account the fact that these planets are 
interacting and are part of a five-planet system. 

3.2. Eccentricity evolution assuming the five-planet system 

3.2.1 . Tides in multiple planets systems 

To compute the tidal interactions we used the tidal force as ex- 
pressed in |Hut| ( [T98T] ) for the constant time lag mod el. We added 
the force in the N-body code Mercury (Chambers] 1999| ). Kaib 
et al. ( 2011| ) showed that the system is expected to be nearly 
coplanar (all the orbits remain within 5° in their simulations). 
Ehrenreich et al.| ( |2012| ) claim that they have observed the sig- 
nature of an extended hydrogen atmosphere of planet b, grazing 
the stellar disk, implying a very low mutual inclination between 
planets e and 'b'. Because planet e is transiting, assuming copla- 
narity implies that all the minimum masses inferred from radial 
velocity are close to the true masses. Kaib et al. (201 1) predicted 
that the plane of the orbits should be misaligned with the star 
due to the influence of the binary companion 55 Cnc B. In this 
study, however, we assume a zero inclination of the system and 
neglect the gravitational influence of 55 Cnc B. Furthermore, 
we consider that the planets have zero obliquity. We consider 
the tidal forces between the star and the planets but we neglect 
the tidal interaction between planets. In this formalism, the total 
tidal force exerted by the star on the i th planet gives an accelera- 
tion for the planet of: 



^M a z,tides — 



W 



- -f (Mfk 2 ,*Rl + Mlk 2 ,tf) 

i 

+ ^[MfcrX°(^*-4-) 

i 

+ Ml<r#}°(Q, - 6i)]eo,i, 



**r,i 



(4) 



The tidal evolution of a single planet around a star is deter- 
mined by the tide raised by the planet on the star - the stellar tide 
- and the tide raised by the star on the planet -the planetary tide. 



The eccentricity damping timescale r ei 
tide is given by r et 



\plan 



due to the planetary 



\plan 



\plan' 



where Rf is the radius of the i th planet, M t its mass, <x; its 
dissipation factor, & 2 ,* its potential Love number of degree 2. r t 
is the distance between the center of star and the i th planet and ft 
its derivative. 0; is the derivative of the i th planet's true anomaly 
0i (the instantaneous orbital angular velocity). The radial vector 
e r j links the center of the star S to the center of the i th planet P/, 

SP 

e r,i ~ ispj • The orthoradial vector egj is such that e r , 



ISP/I ' 



xe ^ = ^> 
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where hi 



M+Mj 
M++Mi 



Y t x Vj is the orbital angular momentum of the 



i th planet, here in heliocentric coordinates. 

Here, we consider that only the two inner planets feel the 
effects of tides. This assumption is justified by the order of mag- 
nitude calculation made in section |37T| Both planets are assumed 
to be in pseudo-synchronous rotation. 

We also add the orbital accel eration due to General Relativity 
( [Kidder|[T995J|Mardling & Lin|[2002| ). Adding this force causes 
a change in the period of the eccentricity oscillations but does 
not change much the equilibrium eccentricity of planet e. 

We performed different simulations with two sets of initial 
conditions. For each set of initial conditions, we varied the plan- 
etary dissipation from 10" 5 x cr p to 100 x cr p . 



The first set of initial conditions is from Dawson & Fabry cky 
( |2010| ) and the second one is from |Endl et al. f2012|). Table |2 



pOTO ) and E ndl et al.| | |20l2) respectively. 



and T able [3] s how the initial c onditions for Dawson & F abrycky 



Planet 



M sin i a (AU) 



u (deg) 



55 Cnc e 


8.26 M e 


0.01560 


0.17 


181 


55 Cnc b 


0.825 Mj 


0.1148 


0.010 


139 


55 Cnc c 


0.171 Mj 


0.2403 


0.005 


252 


55 Cnc f 


0.155 Mj 


0.781 


0.3 


180 


55 Cnc d 


3.82 Mj 


5.74 


0.014 


186 



Table 2. Minimum mass, and orbital elements of the three outer 
planets of the 55 Cnc system.The values of M sin /, a, e and oj 
come from table 10 of Dawson & Fabrycky (2010). 



Planet 



M sin i a (AU) 



co (deg) 



55 Cnc e 


8.37 M e 


0.01544 


0.0 


90 


55 Cnc b 


0.80 Mj 


0.1134 


0.004 


110 


55 Cnc c 


0.165 Mj 


0.237 


0.07 


356 


55 Cnc f 


0.172 Mj 


0.77 


0.32 


139 


55 Cnc d 


3.53 Mj 


5.47 


0.02 


254 



Table 3. Minimum mass, and orbital elements of the five planets 
of the 55 Cnc system.The values of M sin /, a, e and co are from 
Tables 2 and 3 of |Endl et al.|J20l2) . 
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Fig. 7. Evolution of semi-major axis and eccentricity for the 
five planets of 55 Cnc assuming the "nominal" tidal dissipation 
for 55 Cnc e an d the initial conditions come from Dawson & 
Fabrycky ( 2010| ). Top panel: the semi-major axis is plotted in 
thick colored lines, the colored dashed lines represent the dis- 
tance of aphelion and perihelion for each planet. The corota- 
tion radius is indicated by a black dashed line and the thick red 
dashed line is the stellar radius. Bottom panel: the colored full 
lines represent the eccentricities of the five planets, the dashed 
dotted lines represent the secular evolution 55 Cnc e (purple) 
and 55 Cnc b (blue) would experience if alone in the system. 
The thick red dashed line represent the critical eccentricity of 
planet e above which its secondary transit depth is increased by 
more than 10% by tidal heating. 



3.2.2. Results 



Orbits and masses from D awson & Fabrycky] ( |2010| l 

Figure [7] shows the 1 Myr evolution of the 5 -planet system 
starting from the orbits from Tableland for a dissipation factor 
of 1 X cr p for 55 Cnc e. For the first few 10 4 yrs, the eccentricity 
of 55 Cnc e (initially set to 0.17) decreases as expected by the 
secular evolution calculations (dashed-dotted purple line), as if 
it were alone in the system, save a few oscillations. But after 
this initial decay, the eccentricity stops decreasing and oscillates 
around a few 10" 3 , due to the simultaneous influence of tidal 
damping and eccentricity pumping by the other planets. 

Over 1 Myr, the semi-major axis of 55 Cnc e does not vary 
significantly. It decreases from 1.560 x 10" 2 AU to 1.515 x 
10" 2 AU in the first few 10 4 yrs while the eccentricity is still 
high and then decreases to reach 1.512 x 10" 2 AU at 1 Myr. As 
expected from the order of magnitude timescale, the eccentricity 
of 55 Cnc b does not decrease over the integration time. However 
due to the presence of other planets exciting its orbit, 55 Cnc b's 
eccentricity oscillates around 0.016. 



We tested if assuming pseudo- synchronization for 55 Cnc e 
and 55 Cnc b is correct when their eccentricities oscillate around 
their mean value. When 55 Cnc e has the maximum eccentricity 
its pseudo- synchronization period is of 17.949016 hr, and when 
it has the minimum eccentricity, it is of 17.950161 hr. The differ- 
ence between these two values is of < 0.01%, which is negligible 
for dissipation. So the pseudo- synchronization remains a good 
assumption. For 55 Cnc b tidal interaction is weak so it keeps its 
initial rotation period throughout the evolution. We chose here 
to assume pseudo-synchronization. 

To validate our version of the N-body code Mercury with 
tides, we verified it could reproduce the tidal evolution of single 
planet syst ems simulated with th e code used in Bol mont et al. 
(2011) and Bolmo nTet al.| ( |2012| ), which solves the tidal secular 
equations for the semi-major axis, eccentricity (equation[3]of this 
paper) and rotation spin evolution of planet and star. We tested 
the General Relativity effects with Mardling (2007) 's calcula- 
tions for the secular evolution of two-planet systems. However 
as they use the constant phase lag model to compute the tidal 
force, we cannot compare their tidal evolution calculations with 
our simulations, which are based on the constant time lag model. 
The recent literature has debated the applicability of the con- 
stant time lag and constant phase lag models. The constant phase 
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lag model has been shown to be mathematically i nconsistent 
(Efroimsky & Makarov 2013). On the other hand, Greenberg 
( 2009J) showed that the response of a viscoelastic fluid reduces 
exactly to the constant time lag approach in the limit where the 
forcing frequency is smaller than the natural frequency of the 
system and where damping is weak (see his equations 8b and c 
in these limits). A real planet's response is certainly much more 
complex than the one of a viscoelastic fluid, meaning that the 
constant time lag approach must also be used with care. The 
constant time lag model is a linear theory but taking into account 



10.000 
1 .000 k 



3 

< 0.100 

CO 



0.010 r 



0.001 
1 .0000 



0.1000 r 



0.0010 r 



0.0001 



nonlinear terms can result in noticeable changes (Weinberg et al. 
2012). However, we choose this model because we believe that 
it stands on more robust mathematical grounds than the constant 
phase lag model. 

As the dissipation factor of planets is not well known, we de- 
cided to cover a wide dissipation factor range, from 10" 5 x cr p to ° 0.0100 
100 x cr p . As seen in section [51 the dissipation factor affects the 
tidal heating and change the value of the eccentricity required to 
affect the depth of the secondary eclipse. In addition, changing 
the dissipation in 55 Cnc e has two effects on its orbital evolu- 
tion. First, it has an effect on how fast its eccentricity decreases. 
For a dissipation factor of 0.1 x cr p , the eccentricity reaches the 
"equilibrium" in ~ 3 x 10 5 yrs, while for a nominal dissipa- 
tion factor, it does so in ~ 3 x 10 4 yrs, and ~ 3 x 10 3 yrs for 
10 X a p . This is in ag reement with the order of magnitude calcu- 
lation of section [37T] In fact, the eccentricity damping timescale 
r ecc scales as T v so as cr~ l oc (k2, p AT p )~ l , so if the dissipation 
factor is twice the fiducial value r ecc is half the value calculated 
for the nominal dissipation factor. 

Thus for dissipations below 10" 2 x cr p the "equilibrium" ec- 
centricity is not reached in less than 1 Myr because the system 
is still in the circularization period. For dissipation lower than 
10" 2 x <T P , the minimum, average and maximum values of the 
eccentricity of Table [5] are calculated for the whole evolution. 
For dissipation bigger than 10" 1 x cr p , the "equilibrium" eccen- 
tricity is reached in less than 1 Myr. The minimum, average and 
maximum values of the eccentricity are then calculated when the 
eccentricity has reached the "equilibrium". 

The evolution of the system assuming 55 Cnc e has a dissi- 
pation factor of 10" 2 x cr p can be seen in Fig.p] The eccentric- 
ity of 55 Cnc e remains higher than the critical eccentricity for 
~ 6 x 10 5 yrs. For a dissipation of 100 x cr p , the planet 55 Cnc e 
falls on its star in ~ 1.2 x 10 5 yrs. It is therefore likely that the 
planet dissipates less than this value because the probability of 
observing the system in this configuration is low. 

Changing the dissipation factor of 55 Cnc e has also an effect 
on the forcing of its eccentricity by the other planets. The "equi- 
librium" eccentricity oscillates between a minimum and a maxi- 
mum value, which are given in Table[5]for the different planetary 
dissipation factors. 

As expected, when the planetary dissipation is higher the 
forced eccentricity is lower. In Fig.|7j we plot the critical eccen- 
tricity above which the tidal heating has an effect of more than 
10% on the depth of the secondary transit. This value evolves 
with the dissipation. If the planet has a large dissipation fac- 
tor, the eccentricity needed to have a strong tidal dissipation is 
lower. In section |4| we discuss the configurations that produce 
the strongest tidal signature in the transit depth. 
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Fig. 8. As Fig.Mbut assuming a dissipation factor of 10 2 x cr p 
for 55 Cnc e. 
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Fig. 9. As Figure [7J but the initial conditions come from 
|etaL| ( [20T2l ). 



damped to values of a few 10" 3 in less than 10 5 yrs. Table [5] 
also shows that the eccentricity of 55 Cnc e never becomes high 
enough to cause a notable increase in the IR emission of the 
planet. In the absence of significant tidal heating, inefficient re- 
distribution of heat and an albedo lower than 0.55 are required. 



FigfT] shows that the ecce ntricity of 55 Cnc e determined 
by |Dawson & Fabry cky| pOTO}, e = 0.17, is not stable as it is 



Orbits and masses from |Endl et al.^2012) 

Using the orbital distances, masses and exact arguments of 
pericenter proposed by |Endl etaL ( 2012 ) lead to a destabilization 
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10- 5 o> 


10" 4 x o- p 


10" 3 x o- p 


10" 2 x o-p 


10- 1 x o-p 


10° x o-p 


10 1 x o-p 


10 2 x o-p 


min ecc 


1.66 xlO" 1 


1.63 xlO" 1 


1.33 xlO" 1 


1.56 xlO" 1 


7.70 x 10" 5 


5.00 xlO" 6 


6.70 x 10" 5 


1.00 xlO" 6 


mean ecc 


1.70 xlO" 1 


1.68 xlO" 1 


1.52 x 10 1 


6.82 x 10" 2 


1.75xlO" 3 


9.36 x 10" 4 


3.16 xlO" 4 


4.45 x 10" 5 


max ecc 


1.73 xlO" 1 


1.73 xlO" 1 


1.72 x 10 1 


1.71 x 10 1 


3.27 x 10" 3 


1.80xlO" 3 


6.85 x 10" 4 


9.10 xlO" 5 


crit. ecc 


4.21 x 10" 1 


2.71 x 10" 1 


1.31 x 10 1 


4.76 x 10" 2 


1.52 xlO" 2 


4.82 x 10" 3 


1.52xl0" 3 


4.82 x 10" 4 



Table 4. Minimum, maximum and mean value of the eccentricity of 55 Cnc e obtained in the simulation of the system proposed by 
Dawson & Fab rycky (2010) for different dissipation factor for 55 Cnc e. The critical eccentricity needed to affect the planet is also 
indicated for each case. 





lO-V, 


10" 4 x o-p 


10" 3 x o-p 


10" 2 x o-p 


10- 1 x o-p 


10° x o-p 


10 1 x o-p 


10 2 x o-p 


min ecc 


1.07 xlO" 4 


1.22 xlO" 4 


1.05 xlO" 4 


3.10 xlO" 5 


2.50 x 10" 5 


7.16 xlO" 4 


6.00 x 10" 6 


0.00 


mean ecc 


7.55 x 10" 3 


5.78 x 10- 3 


4.36 x 10" 3 


4.33 x 10" 3 


3.27 x 10" 3 


4.98 x 10" 3 


8.16 xlO" 4 


1.74 xlO" 5 


max ecc 


1.30 xlO" 2 


1.11 xlO" 2 


1.03 xlO" 2 


1.01 x 10" 2 


9.38 x 10" 3 


1.04 x 10" 2 


1.52 x 10- 3 


4.54 x 10" 4 


crit. ecc 


4.21 x 10" 1 


2.71 x 10" 1 


1.31 x 10" 1 


4.76 x 10" 2 


1.52 xlO" 2 


4.82 x 10" 3 


1.52 x 10" 3 


4.82 x 10" 4 



Table 5. Minimum, maximum and mean value of the eccentricity of 55 Cnc e obtained in the simulation of the system proposed by 
Endl et al. ( 2012| ) for different dissipation factor for 55 Cnc e. The critical eccentricity needed to affect the planet is also indicated 
for each case. 



of the system. We simulated the system with the N-body code 
Mercury without the tidal forces and without General Relativity 
and the outcome is the ejection of a planet. We also simulated 
the system without tides but with General relativity and the out- 
come is also the ejection of a planet. However, if we simulate 
the sys tem with the orbital distances and masses from |Endl et al. 



(2012) but using random arguments of pericenter (knowing that 



the values they provide are affected by significant uncertainties), 
we usually find that the system is stable on 10 7 yrs timescales. 

Fig. [9] s hows the simulati on using orbital elements and 
masses from Endl et al. (2012). The important differences be- 
tween [EndTetlipl^^ for the 
dynamical evolution of the system mainly li e in the value s of the 
eccentricities of b and c. In Dawson & Fabrycky ( |2010 ), b and 
c have small eccentricities: 0.01 for b and 0.005 for c, while in 



End l et al.| p012) b has a small eccentricity: 0. 004 but c has a 
high eccentricity: 0.07. In |Dawson & Fab rycky (2010) after the 
damping of its initially high value, the eccentricity of e depends 
mainly on the eccentricity o f b and in a weaker way on the ec- 
centricity of c. However for Endl et al. ( 2012| ), we observe that 
c excites b from the initial value of 0.004 to the mean value of 
-0.02. So e is excited by a more e ccentric planet b than was the 
case for pawson & F abrycky ( 2010| ), and in a weaker way by a 
more eccentric planet c. 

The eccentricity of 55 Cnc e cannot be zero because of the 
presence of the other planets, it oscillates around 4.98x 10" 3 . The 
eccentricity of b also increases to oscillate around 2.20 x 10" 2 . 
However, during the evolution of the five planets, the system 
frequently passes within the error bars of the current configura- 
tion derived by (Endl et al. 2012| ). For these configurations, the 
eccentricity of 55 Cnc e is higher than the critical eccentricity 
for tidal heating. This simulation therefore provides dynamical 
evolution snapshots that are simultaneously consistent with the 
system properties derived from radial velocities, the secondary 
eclipse depth observed with Spitzer and an observable tidal heat- 
ing. 

As seen in Table [5] and Fig. [9] the mean eccentricity of 
55 Cnc e is also higher than the critical value. For 10 x cr p , the 
maximum value of the eccentricity is above the critical eccen- 
tricity while for 0.1 x cr p the eccentricity always remain below 
the critical value. 

In Fig [3] (airless case) we can see that the upper limit on 
the albedo is 0.75 for the mean eccentricity from this simula- 



tion. For the maximum eccentricity, this upper limit is as high as 
0.85. Interestingly, these high values of the albedo match those 
inferred for the exoplanet Kepler 10 b. 

In the isothermal case (Fig [4]) the upper limit on the albedo 
is 0.1 for the mean eccentricity from this simulation. For the 
maximum eccentricity, this upper limit is 0.7. No solutions are 
found for 0.1 x cr p and a marginal match is found with 10 x cr p 
withA<0.05. 



4. Tidal dissipation, transit depth and eccentricity 

Fig. [10] sums up what has been discussed so far. It shows the 
range of simulated eccentricities of 55 Cnc e for an evolution 
of 1 Myr superimposed with the albedo-dependent region in the 
eccentricity-planetary dissipation parameters space correspond- 
ing to the two following conditions: 



^p,tides/^p,tot > 0.1 

103 ppm < Fp.tot/F* < 159 



ppm 



(5) 



The first condition is that the tidal flux is more than 10% the 
total flux of the planet. The second condition corresponds to the 
observability constrain of a transit depth of 13 1 ± 28 ppm, where 
F+ is the flux of the star calculated - as the flux of the planet - 
at a distance of 12.34 pc (D12) in the Spitzer band. 

The configurations for which the two conditions are met are 
therefore where the thick colored curves representing the output 
of our Mercury code intersect the shaded areas. For pawson &| 
Fabrycky] ( |2010| ), the overlap occurs from a dissipation of a few 
~ 10~* xcr p to a dissipation of a few 10" 2 x cr p . However, as the 
eccentricity of the planet has been constrained by observations 
(D12) to be smaller than 0.06, the lower limit of the dissipation 
we can infer is about 8 x 10" 3 cr p - where the dashed dotted 
orange line intersects the horizontal black dashed line. For dissi- 
pations lower than 10" 2 xcr p , the eccentricity of 55 Cnc e has not 
reached the "equilibrium" eccentricity, and the planet is still in 
its circularization period after 1 Myr. The eccentricity decreases 
slowly with small oscillations as seen in Fig. [8] So for dissipa- 
tions lower than 10" 2 x cr p , the minimum eccentricity plotted in 
Fig.flOHs the eccentricity of 55 Cnc e at t = 1 Myr and the max- 



imum eccentricity is the high initial one 0.17. It means that for 
these configurations, the two conditions of equation [5] are met 
for a short period of time only. For example, for a dissipation 
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Fig. 10. Eccentricity vs dissipation. The thick colored lines cor- 
responds to the minimum, average and maximum eccentricity 
- bottom full line, middle dashed line and top full line - of 
55 Cnc e obtained wit h the simulations of Dawson & F abrycky 
( |2010| ) in purple, and |Endl et al.| ( [20T2| ) in blue. The horizontal 
black dashed line corresponds to the upper limit from observa- 
tions (e = 0.06, D12). The color bands correspond to the region 
where the two conditions of equation [5] are met: where the tides 
contribute to more than 10% of the planet's flux and where the 
resulting transit depth is between the observational constraints. 
The light orange band corresponds to an albedo of and the red 
one to an albedo of 1 . The orange dashed dotted lines delimit the 
band corresponding to an albedo of 0.65. 



5. Discussion 

This work shows that the albedo cannot be constrained directly 
from the eclipse depth due to the influence of tidal heating 
on the radiative budget. For a given dissipation constant there 
is, therefore, a degeneracy between albedo, eccentricity and 
the assumed heat redistribution scenario. This degeneracy can 
partially be broken by using precise orbital photometry, as 
shown in fig. [T] A measure of the eclipse depth and/or the 
phase curve a different wavelengths would provide an additional 
constraint on the temperature map (as shown for instance by 
Maurin et al., 2012) and thus on the heat redistribution as well 
as on the relative contributions of tidal and radiative heating. 
A measure of the secondary eclipse and/or the phase curve at 
visible wavelengths that are not affected by thermal emission 
would provide a direct measurement of the albedo and would 
also help breaking this degeneracy. The ratio between the day- 
side emission of the planet and the star is however very small: 
Ax 17.5 ppm and should be measured over a broad wavelength 
range in order to infer the Bo nd albedo. The p hotometric 
survey of 55 Cnc with MOST ( |Winn et all |2011| ) reveals a 
modulation of the luminosity of the system (in the 350 - 700 nm 
window), at the orbital period of the planet and whose phase is 
consistent with a phase modulation. However, the amplitude of 
this modulation, -170 ppm is nearly 10 times larger than what 
could be attributed to the planet. In addition the light curve does 
not exhibit a significant drop when the secondary eclipse o ccurs. 
This points to a modulation of the stellar luminosity itself. Winn 
et al.| ( |2011| ) suggest that that part of the stellar hemisphere 



facing the planet could be slightly fainter than the rest of the star 
due to magnetic interaction between the planet and star. Such a 
stellar modulation would also make extremely challenging the 
extraction of the planetary phase modulation. How better is the 
situation in the infrared still has to be evaluated. 



of 10" 2 x cr p , the two conditions are met for - 6 x 10 5 yrs after 
the proposed present configuration (Dawson & Fabrycky, 2010). 
We do not know the previous state of the system, but it is likely 
that with such low dissipation the ancient eccentricity was sig- 
nificantly higher than its present value. 

The most favorable configuration is for a planet with A ~ 
0.65 as can be seen in Fig. [5] This albedo corresponds to the 
band delimited by orange dashed dotted lines in Fig. 10 Fig. [3] 
shows that if the tides contribute to more than 10% to the planet's 
flux, and the eccentricity of the planet is less than 0.06 and the 
transit depth is in agreement with the observations, the albedo of 
the planet would be between ~0.5 and 0.7. 



For |Endl et al.| ( |2012] ), a consistency between the two condi- 
tions and the eccentricity range obtained with Mercury simula- 
tions is found for dissipations around 1 x cr p . Depending on the 
dissipation, the two conditions are in agreement with the whole 
range of albedo. Contrary to the case of |Dawson & Fabrycky| 
( 2010| ), the eccentricity oscillates around the "equilibrium" ec- 
centricity for all the dissipation range values. So over time the 
system will be periodically in configurations for which the two 
conditions are met. Over an oscillation of the eccentricity, the 
two conditions are met when the eccentricity is the higher. For 
example, as can be seen in Fig. 10 for a dissipation of 10 x cr p , 
the only configuration which would be in agreement with the 
two conditions is when 55 Cnc e reaches the maximum value of 
the eccentricity oscillations, and it would mean that the planet 
has an albedo of 0.65. 



Our calculations assume that the internal heat flux is released 
uniformly over the planet. The existence of one or a few hot 
spots could significantly change the results. As an extreme case, 
if all the heat is released in the dark hemisphere, it has no effect 
on the dayside emission and thus a weak effect on the secondary 
eclipse depth. This would however have a strong signature on 
the orbital photometry. 

In this study, we modeled the orbital photometric signature 
of planet e as if it was alone in the system. |Kane et aL] ( |2011| ) 
showed that, at visible wavelengths, the modulation due to 
planet 'b' may not be negligible compared with that of planet 
e, depending of course on the respective albedos of planets e 
and 'b'. Characterizing the properties of planet e by orbital 
photometry thus requires the ability to distinguish between the 
two components, which is an additional challenge. At thermal 
wavelengths, the situation should be different because planet 'b' 
is a Jovian planet and atmospheric circulation tends to attenuate 
the day/night temperature contrast and thus the amplitude of 
thermal phase curve. To estimate this temperature contrast we 
can compare radiative and dynamical timescales. The radiative 
cooling timescale is given by tr = ^^ where cp is the specific 
heat of H2, P is the pressure at the considered level and T is the 
bolometric brightness temperature. Assuming a superrotating 
atmosphere (as that of Venus or as predicted for strongly 
irradiated atmospheres (see for instance Showman et al., 2009 
and Selsis et al., |2011| )), a relevant dynamical timescale is given 
by r D = R/v where R is the planetary radius and v is the average 
of the zonal wind. Let us use P = 0.1 bar (a typical photosphere 
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for infrared radiation), T=700 K (the equilibrium temperature 
of the planet for A = 0), R = Rj up . Using 3D hydrocodes, ISelsis 



et al.| ( |2011| ) found v ~ 200m/s for the mean zonal wind speed 



for a 10 bar CO2 atmosphere and an equilibrium temperature 
of 390 K, while |Showman eFaT] ( [2009] ) found one order of 
magnitude faster winds for HD 189733b (T eq = 1200 K). 
Assuming 200 < v < 2000 m/s yields 1 < t r /t d < 10, which 
implies an efficient heat redistribution between the day and 
night hemisphere and thus a rather flat phase curve. The change 
of equilibrium temperature between periastron and apoastron is 
less than 5 K, a difference that should be completely averaged 
out by the inertia of the atmosphere. The contribution of planet 
'b' to the photometric variations of the system should thus be 
negligible at thermal wavelengths. 

Observations in several infrared bands will be necessary to 
establish the radiative budget of the planet and demonstrate the 
existence of an excess of luminosity due to tidal heating. In par- 
ticular, in the presence of an atmosphere, the measured bright- 
ness temperature can correspond to the physical temperature at 
any atmospheric level. This temperature can locally exceed the 
highest temperature that can be found on an airless planet with 
a null albedo. This can happen in two cases. First, if there is 
a strong greenhouse warming and if the observed band corre- 
sponds to an atmospheric window that probes the surface of the 
planet. This situati on can be illustra ted by the thermal phase 
curves modeled by Selsisetal. ( 2011| ) for C02-rich atmospheres. 
In their Fig. 6, for instance, one can see that the 8.7 jim band 
(an atmospheric window) exhibits fluxes that are higher in the 
presence of an atmosphere. However, expected components of 
a planetary atmosphere, like CO2, CO, SO2, H2S, CS2, have 
strong absorption bands in the IRAC2 spectral window, mak- 
ing unlikely for this band to probe the deep and hot atmospheric 
layers heated by greenhouse effect. An "excess" of brightness 
temperatures can also be observed in an opaque band probing 
high stratospheric layers subjected to a strong thermal inversion. 
In the present case this would imply that the atmosphere is suf- 
ficiently opaque over the whole 4 to 5 yum for the emission to 
come from very low pressure layers. Indeed, on the day side of 
a synchronously rotating planet, the convective plume powered 
by the stellar irradiation is drafted to very high altitudes and an 
inversion appears possible only in the upper atmosphere. 

In a recent work, Maurin et al. (2012) proposed to use multi- 
wavelengths orbital photometry to derive the albedo, radius and 
inclination of rocky planets on circular non-transiting short- 
period orbits. The present study shows that applying this tech- 
nique requires to insure that tidal heating is not affecting the ra- 
diative budget of the planet. Planets subjected to the influence 
of a massive and eccentric companion could have a small non- 
measurable eccentricity sufficient to produce a significant tidal 
heating. 

In our simulations, we assume the pseudo- synchronization 
of the planet spin. For rocky planets, the very existence of 
such state has been recently contested by Ma karov & Efroimsky| 
( 2012| ), who suggest that the only likely rotation configurations 
are spin-orbit resonances. Taking into account a more realistic 
rheology for the rocky planets, they show that the seemingly 
lowest-energy state of pseudo-synchronization is actually unsta- 
ble. In the present study, the eccentricities remain small enough 
so that the difference between the 1:1 spin-orbit resonance and 
pseudo-synchronization has a negligible effect on dissipation, 
at least after a few 10 4 yrs. The evolution during the first few 
10 4 yrs of the simulation of the Dawson & Fabry cky ( 2010| 
system would, however, be different. Using the constant time 



lag model, Rodriguez et al. ( 2012| ) showed that even though 
spin-orbit resonances are possible for large eccentricities, plan- 
ets could be in pseudo- synchronization for eccentricities smaller 
than 0.1. This study takes into account the oblateness of the 
planet in a most accurate way than in our model but does not 
in corporate a realistic rheology f or the terrestrial planets such as 
in[ 



Makarov & Efroimsky (2012). 



The initial conditions of our simulations are the present val- 
ues of semi-major axis and eccentricity of the planets. However, 
planet e must have formed farther out and migrated in due to tidal 
evolution. Tides are a dissipative process so it is delicate to inte- 
grate the evolution back in time, the solution diverges quickly to 
a huge semi-major axis and an eccentricity close to 1. To inves- 
tigate a possible formation distance from the star would require 
to explore a range of semi-major axis initial conditions to try 
reproduce present values. Our simulations are time demanding, 
so exploring a huge parameters space - including the planetary 
dissipation- in order to find a appropriate configuration would 
require a lot of time. So, we chose here not to treat past evolu- 
tion, but investigate the present dynamics of the system. 

In this paper, we show the influence of the eccentricity of the 
outer planets on the eccentricity of the inner planet. In the case 
of End l et al. ( 2012| ), the presence of the eccentric planet 'c' can 
increase the eccentricity of e from the supposed observational 
value of to 5 x 10" 3 . The presence of the other planets also 
influences the eccentricity of the inner planet, so the values of 
the inner planet equilibrium eccentricity could change if more 
planets exist in the system as might be the case according to 
|Raymond etaT]{|2008). 



6. Conclusions 

Our study shows that the orbit of exoplanet 55 Cancri e cannot 
be fully circularized due to the influence of the other planets in 
the system. The forced eccentricity depends on the unknown 
dissipation constant and on the detailed set of orbital elements 
for the system. Using an Earth- like dissipation and orbital 
parameters from En dl et al.| ( |2012| ), we find an eccentricity that 
oscillate around a mean value 5 x 10" 3 , up to 8.5 x 10" 3 . The 
resulting tidal heating is sufficient to affect the thermal emission 
that has been observed at secondary eclipse with Spitzer/IRAC2 
by |Demory et al.|(2012| ). 



With e = (and thus no internal heating) we find the same 
constraints on the planetary albedo as |Demory et al. ( 2012| ): 
A < 0.55 in the absence of heat redistribution and no match in 
the case of a uniform temperature. Including tidal heating the 
upper limit on the albedo becomes 0.9 and 0.7 without heat 
redistribution and in the uniform case, respectively. With these 
new constraints, both 55 Cnc e and Kepler 10 b can have a high 
albedo (A > 0.5). 

Assuming an Earth-like dissipation allows us to use the 
Spitzer measurement to put an upper limit on the eccentricity: 
e < 0.015, which is significantly lower than the observational 
constraint from D12 (e < 0.06). The actual dissipation is 
unknown but future observations (in particular muti- wavelength 
orbital photometry) will have the ability to better constraint 
the contribution of tidal dissipation and thus the eccentricity. 
Alternatively, a more accurate measurement of the eccentricity 
and of the inclination of the system would allow us to constrain 
the dissipation. 
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This shows the importance tidal heating on the radiative 
budget of close-in exoplanets, in particular in compact systems 
that include eccentric planets. This case illustrates how the dy- 
namics of a whole planetary system can influence the retrieval 
of the intrinsic properties of one particular planet (here, for 
instance, the surface albedo). 
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